Abstract We discuss the capability of a third-generation ground-based detec- 
tor such as the Einstein Telescope (ET) to enhance our astrophysical knowl- 
edge through detections of gravitational waves emitted by binaries including 
intermediate-mass and massive black holes. The design target for such instru- 
ments calls for improved sensitivity at low frequencies, specifically in the <~ 1- 
10Hz range. This will allow the detection of gravitational waves generated in 
binary systems containing black holes of intermediate mass, ~ 100 lOOOOM©. 
We primarily discuss two different source types — mergers between two in- 
termediate mass black holes (IMBHs) of comparable mass, and intermediate- 
mass-ratio inspirals (IMRIs) of smaller compact objects with mass ~ 1-1OM 
into IMBHs. IMBHs may form via two channels: (i) in dark matter halos at 
high redshift through direct collapse or the collapse of very massive metal-poor 
Population III stars, or (ii) via runaway stellar collisions in globular clusters. 
In this paper, we will discuss both formation channels, and both classes of 
merger in each case. We review existing rate estimates where these exist in 
the literature, and provide some new calculations for the approximate num- 
bers of events that will be seen by a detector like the Einstein Telescope. These 
results indicate that the ET may see a few to a few thousand comparable-mass 
IMBH mergers and as many as several hundred IMRI events per year. These 
observations will significantly enhance our understanding of galactic black-hole 
growth, of the existence and properties of IMBHs and of the astrophysics of 
globular clusters. We finish our review with a discussion of some more specu- 
lative sources of gravitational waves for the ET, including hypermassive white 
dwarfs and eccentric stellar-mass compact-object binaries. 
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1 Introduction 

The Einstein Telescope (ET), a proposed third-generation ground-based gra- 
vitational-wave (GW) detector discussed in greater detail elsewhere in this 
volume, will be able to probe GWs in a frequency range reaching down to 
~ 1 Hz 51,36]. This is lower than the limit of ^40 Hz available to current 
ground-based interferometric GW detectors such as LIGO, Virgo, and GEO- 
600 or the ~ 10 Hz limit that could be reached by their second generation 
|105l [3ll4"3"] . On the other hand, GWs in the range above ~ 0.1 Hz will not be 
accessible to the planned Laser Interferometer Space Antenna (LISA, [IS]), 
which will have sensitivity in the ~ O.lmHz-O.lHz range. The frequency range 
determines the typical masses of coalescing binaries that could be detected 
by an interferometer; for example, the frequency of GWs emitted from the 
innermost stable circular orbit of a test particle around a Schwarzschild black 
hole of mass M is « 4400 Hz (M Q /M). The Einstein Telescope will therefore 
probe sources with total masses of hundreds or a few thousand M Q which 
are less likely to be detected by LISA or the current ground-based detectors. 
This places the ET in a position to make complementary observations to LISA 
and LIGO/Virgo/GEO-600 and to carry out unique searches for several very 
exciting source types, particularly those involving light seeds of massive black 
holes and intermediate-mass black holes. 
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There is a significant body of evidence that massive black holes (MBHs) are 
generically found in the centers of massive galaxies [31]. These MBHs merge 
during mergers of their host galaxies, and such mergers therefore trace the 
history of structure formation in the universe. Gravitational waves emitted 
during the mergers of MBHs with masses in the ~ 5x 10 4 Mq-5 x 1O 7 M will 
be detectable by LISA; dozens of detections could be made during the LISA 
mission |101j . According to some predictions, these massive black holes grow 
from light seeds of ~ 100 M through accretion and mergers [55111 121H02"] . 
The typical frequencies of gravitational radiation emitted during the mergers 
of such systems will fall in the 0.1-10 Hz range, however, and will only be 
accessible to GW detectors sensitive in that range. The Einstein Telescope 
may be able to detect tens of such sources, determining their masses to an 
accuracy of a few percent and the luminosity distances to < 30% |100ll3T)] . 

Additionally, globular clusters may host intermediate-mass black holes 
(IMBHs) with masses in the ~ 100 - 1000 M Q range (see [75UT2] for re- 
views). Intermediate-mass-ratio inspirals (IMRIs) of neutron stars or stellar- 
mass black holes into these IMBHs could be detected by the second generation 
of ground-based detectors [70]; however, the Einstein Telescope should be able 
to detect far greater numbers of events, up to as many as several hundred per 
year, and these events will range to higher IMBH masses. If the binary fraction 
in a globular cluster is sufficiently high or two globular clusters hosting IMBHs 
merge, an IMBH-IMBH binary can form and then coalesce, emitting gravita- 
tional waves [35J . The Einstein Telescope could detect as many as thousands of 
such events, although, given the present uncertainty about the very existence 
of IMBHs, all such estimates must be viewed with a great deal of caution. 

The Einstein Telescope may also be able to detect a number of other, more 
speculative, sources. These include the inspirals of stellar-mass black holes into 
IMBHs that may reside at the centers of dwarf galaxies, although we do not 
expect a significant rate of detectable signals of this type. There is also the 
possibility that ET will detect orbiting white dwarfs near the upper end of 
their allowed mass range, or eccentric compact object binaries. 

This paper is organized as follows. In Section [2] we discuss the method- 
ology for the event-rate calculations. We describe the adopted detector and 
network models, the formalism for estimating the signal-to-noise ratio and the 
waveform families used in the analysis. We then discuss in detail several types 
of GW sources of particular relevance to the ET. In Section [3] we consider 
sources involving IMBHs in globular clusters, both intermediate-mass-ratio 
inspirals and IMBH-IMBH coalescences. In Section [4] we focus on light mas- 
sive black holes. We discuss how mergers between galaxies at high redshift 
can lead to IMBH-IMBH binaries detectable by ET, provided the black hole 
seeds in the galaxies are light. We also describe how light massive black holes 
in the centres of dwarf galaxies could act as IMRI sources. In Section [5j we 
discuss several speculative sources, including hypermassive white dwarfs and 
eccentric binaries. We finish with a discussion, in Section [6j of some of the 
potential scientific implications of ET observations of these sources. Section [JJ 
provides a brief summary. 
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2 Methodology for event-rate calculations 

2.1 The Einstein Telescope configuration 

The design target for the Einstein telescope is a 10km scale interferometer, 
with a factor of ~ 10 increase in sensitivity over Advanced LIGO, and im- 
proved sensitivity at low frequencies. The ET design also calls for the abil- 
ity to measure polarisation at a single site, which requires at least two non- 
coaligned coplanar detectors at the site. The currently favoured configuration 
is a triangular facility, with 10km long arms, and containing three independent 
detectors with 60° opening angles, as this has lower infrastructure costs and 
slightly better sensitivity than two right angle detectors placed at 45° to one 
another [35] . We refer to this triangular design as a "single ET" . In Figure [I] 
we show the target ET noise curve, labelled "ET baseline" [5T]. This noise 
curve is for a single right-angle interferometer with the ET design sensitiv- 
ity. Unless otherwise stated, signal-to-noise ratios (SNRs) etc. will be quoted 
for this configuration. The sensitivities of one 60° interferometer, two right- 
angle interferometers and a single ET are changed relative to this by factors 
of v3/2, y/2 and 3/2 respectively The "ET baseline" design has recently been 
superseded by the curve labelled "ET B" in Figure [T] but we have checked 
that this change does not significantly affect our results, since the noise curves 
are largely similar in the frequency range where massive systems accumulate 
most of their SNR. Figure [l] also shows an alternative 'xylophone' configu- 
ration for ET that was described in [55]. Such a noise curve is realised by 
operating two detectors within the same vacuum system, one optimised for 
low-frequency sensitivity and the second for high frequency sensitivity. The 
net effect on the composite noise curve is an improved sensitivity near 10Hz. 
The ET design and noise curve has not yet been finalised, nor what fraction 
of time ET would spend in high-frequency and low-frequency operation under 
the xylophone configuration. We will see in Figure [4] in Section |4.1| that the 
xylophone mode is to be preferred for the detection of black hole binaries in 
the 100-1000M Q range. 

In Section |4.1| we will quote results for the parameter estimation accuracies 
that are achievable for mergers of light seeds of massive black holes detected by 
ET. As the events are short lived, parameter estimation requires the existence 
of a network containing multiple detectors. The results we will present are 
taken from [100| l39 j . in which four third-generation network configurations 
were considered — (i) one ET at the geographic location of Virgo, plus a 
second right-angle 10km detector at the location of LIGO Hanford or Perth 
(Australia); (ii) as configuration (i) plus a third 10km detector at the location 
of LIGO Livingston; (iii) as configuration (i) but with the Hanford/Perth 10km 
detector replaced by a second ET; and (iv) three ETs, one at each of the three 
sites. 
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2.2 Signal-to-noise ratio 



In Sections |3.1| and |3.2| we will quote signal-to- noise ratios for events detected 
with ET. In this section, we describe how these were calculated. The signal- 
to-noise ratio (SNR) p for a waveform h(t) measured by a single detector with 
one-sided noise power spectral density S n (\f\) is given by p 2 = (h\h). Here 
(a | b) is the noise- weighted inner product 

where d(/) is the Fourier transform of the waveform a(t) and * denotes the 
complex conjugate. 

We define the horizon distance -Dhor as the distance at which an optimally 
oriented, overhead source produces an SNR of 8. The actual gravitational- 
wave emission is not isotropic and nor is the response of the detector. We 
can define the average range as the radius of a sphere whose volume is equal 
to the true (non-spherical) volume in which inspiral sources can be detected 
with an SNR greater than 8. For uniformly distributed, randomly oriented 
sources, and without applying corrections for cosmological redshift, this range 
is related to the horizon distance by D = Dhor/2.26 [33]. In general, the 
relationship between D and Dhor for sources at cosmological distances will 
depend on the cosmology, the redshift dependence of the source distribution 
and even on the intrinsic luminosity of the sources. In this paper, we will ignore 
these issues and use the usual 1/2.26 correction factor, since a more rigorous 
calculation is not currently available. We expect that the error that arises from 
this approximation will be of the same order as the SNR uncertainties that 
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arise from inaccuracies in the waveform model, which we will discuss in the 
next section, but emphasise that this must be verified by a proper calculation 
in the future. 

For multiple detectors, the inner product in equation should be replaced 
by sums of inner products over individual detectors. Thus, the existence of N 
identical interferometers will increase the range by a factor of yN over the 
single-interferometer range if the network SNR threshold is fixed. 



2.3 Waveforms 

In this paper, we will focus on two different types of source — binaries consist- 
ing of two intermediate mass black holes of comparable mass; and intermediate- 
mass-ratio inspirals (IMRIs) of stellar mass compact objects (neutron stars or 
black holes) into IMBHs. We model the gravitational waves emitted by these 
systems in different ways. For comparable mass binary systems at the upper 
end of the detectable mass range, a significant amount of energy is radiated 
during the merger and ringdown phases and so it is important to include these 
in waveform models for signal-to-noise ratio calculation and parameter estima- 
tion. The recent advances in numerical relativity have allowed the construction 
of hybrid waveform models that include inspiral, merger and ringdown in a 
self-consistent way in a single template. There are two models currently avail- 
able, both of which are for systems containing non-spinning black holes. The 
assumption of zero spin and the intrinsic waveform uncertainties lead to un- 
certainties in the SNRs at the level of a few tens of percent, which we will 
discuss later. However, the corresponding uncertainties in the event rates are 
much smaller than the typical uncertainties that arise from the astrophysics, 
since the lower limit on the number of intermediate-mass black holes in the 
Universe is zero. 

The non-spinning "phenomenological" inspiral, merger, ringdown model 
(NSphenom) 013], is constructed by taking a simple ansatz for the waveform 
model, inspired by post-Newtonian theory, and fitting the coefficients in this 
model to numerical relativity simulations. The phenomenological model has 
recently been extended to spinning black holes with non-precessing spins [6]. 
The "effective-one-body, numerical relativity" (EOBNR) model [23] uses post- 
Newtonian expressions to model the inspiral radiation, which are matched 
onto fits to numerical relativity simulations for the merger radiation and then 
onto analytic expressions for the quasinormal mode ringdown radiation. The 
EOBNR waveforms are constructed to match perturbative results in the test- 
particle limit where the mass-ratio tends to zero. Denoting the mass of the most 
massive object in the binary by mj and the less massive object by 7712, the 
NSphenom and EOBNR models give waveforms as a function of the total mass 
of the binary, M = m\ + 7712, the reduced mass ratio, 77 = m\m2/M 2 and the 
time of merger, to. Including the detector response introduces six additional 
extrinsic parameters expressing the relative location and orientation of the 
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source and detector: the distance, the two sky-location angles, two binary 
orientation angles, and the phase at some fiducial time, e.g., to. 

In Section |4.1| we quote results for the SNRs of comparable mass IMBH 
mergers detected by ET. These SNRs were computed in [100] using the NSphe- 
nom model. ET SNRs for the same systems were recomputed in [39 using 
both models and the results were found to agree to ~ 20%, with the NSphc- 
nom predictions being higher for low mass ratios, r\ ~ 0.16, and the EOBNR 
predictions being higher for high mass ratios, r\ rj 0.25. The two waveform 
families are constructed in different ways and have been matched to different 
numerical relativity simulations, so it is not surprising that the results they 
predict differ. The level of difference provides a guide to how much the SNRs 
computed using either waveform family will differ from the true SNRs of these 
systems that nature provides. This ~ 20% SNR uncertainty from the waveform 
model must be compared to other uncertainties, such as the omission of spin 
from the waveform model and uncertainties in the intrinsic astrophysical rates. 
The effect of spin on the average SNR will also be at the ~ 10% level |69| . 
although it is more significant for the highest mass systems as it can bring oth- 
erwise undetectable sources into the frequency band of the detector. However, 
the net effect will be to increase the rate of detectable events so using non- 
spinning models can be thought of as conservative in that regard. A 50% SNR 
uncertainty leads to an uncertainty in the event rate of a factor of ~ 3, but 
typical uncertainties in the astrophysical rates are an order of magnitude. The 
SNR could be out by as much as a factor of 2 before it would be comparable to 
these astrophysical uncertainties, so we consider the waveform uncertainties as 
fairly negligible. Recent comparisons with additional numerical data suggest 
that the errors in the EOBNR waveforms are significantly smaller than those 
in the NSphenom waveforms, at least for equal-mass sources [35]. For this rea- 
son, we use the EOBNR model to compute the new results presented in this 
paper, specifically the SNRs for IMBH-IMBH mergers in globular clusters, 
which we will describe in Section 13.21 

The NSphenom and EOBNR waveforms have been matched to numerical 
relativity simulations, but only for mass ratios of 1:4 and higher. Computa- 
tional requirements suggest it is unlikely that numerical simulations using cur- 
rent techniques will go beyond mass ratios of ~1:10 in the near future, although 
this may be possible using innovative new approaches. Post-Newtonian theory 
also breaks down once the mass-ratio becomes too extreme. For very extreme 
mass ratios, r\ ~ 10 _6 -10 -4 , gravitational waveforms can be computed using 
black hole perturbation theory [89 , in which the smaller object is regarded 
as a perturbing field of the background spacetime of the larger object and 
radiation reaction is described in terms of the 'self-force'. Significant progress 
has been made over the past few years in self-force calculations, which has 
led to the calculation of the self-force for circular orbits in the Schwarzschild 
spacetime [IT] , including the shift in the location of the innermost-stable- 
circular-orbit (ISCO) that results from the action of this force [12j. However, 
even at mass ratios of ~ 10~ 5 , the terms that are missing in the first-order 
self-force formalism are estimated to have a marginal effect on the phasing of 
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waveforms for LISA sources [56] . The mass ratios for typical IMRI sources for 
ET lie somewhere between these extremes, being typically ~ 0.001-0.1. In this 
regime neither post-Newtonian nor pcrturbative waveforms will be adequate 
on their own to model the true waveforms [7T]. More research is needed to 
devise waveforms that are suitable for bridging this gap. In the meantime we 
must make do with the available waveforms, with the understanding that these 
are not completely accurate. 

In Section [XT] we will present ET event rates for IMRI sources. These were 
estimated using SNRs computed with EOBNR waveforms, since these should 
be accurate both for rj 0.25 and in the limit 77 — > 0. To check the validity of 
the results, we also computed SNRs using the NSphenom waveforms and using 
perturbative waveforms for circular and equatorial inspirals, as described and 
tabulated in [33], and as used to estimate SNRs for LISA in [35]. The latter 
waveform model includes only the inspiral phase, and so we compared those re- 
sults to the inspiral contribution to the EOBNR and NSphenom SNRs. For the 
early portion of the inspiral, significantly before the last stable orbit, we found 
that the SNRs obtained with the three waveform families agreed to within 
~ 10%; this difference is of a comparable magnitude to the effect of omitting 
relativistic corrections in post-Newtonian inspirals. We find, however, that for 
more massive systems, when merger and ringdown contribute a significant 
portion of the SNR, estimates from the NSphenom and EOBNR waveforms 
differ very significantly, with NSphenom waveforms predicting SNRs that are 
greater by more than an order of magnitude. In Figure [2] we show how the 
SNRs computed using the two models vary as a function of 77 for a fixed lumi- 
nosity distance and a fixed redshifted total mass of M z — 500M Q . This figure 
clearly shows the large difference between the two models for small 77. Our 
theoretical expectation is that the energy emitted during the ringdown should 
scale as rj 2 x M, so that the ringdown SNR should scale as 77 for a fixed total 
mass. However, the SNR predicted by the ringdown portion of the NSphenom 
waveform scales roughly as y/rj, and is therefore significantly over-predicted 
for small 77. For this reason we used EOBNR waveforms, which exhibit the 
correct scaling, for estimating IMRI SNRs, with the understanding that there 
is a clear need for more careful and accurate modelling of IMRI radiation 
in the future. This is essential not least because detection of these systems 
will almost certainly rely on matched filtering, for which accurate template 
waveforms are needed. 



3 IMBH sources in Globular Clusters 

A particularly exciting possibility for GW astronomy is the observations of 
intermediate-mass black holes (IMBHs) in globular clusters. Unlike stellar- 
mass and supermassive black holes, IMBHs have not been identified dynam- 
ically, that is, through the observation of one or more orbiting objects and 
the application of Kepler's law. Thus, searches for their existence have re- 
lied on circumstantial evidence (for reviews and recent results see [75 |ll09[ l83 ] 
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Fig. 2 Signal-to-noise ratio as a function of mass ratio, r), computed using the NSphenom 
and EOBNR waveform families. SNRs are quoted for systems at a luminosity distance of 
6.61Gpc (corresponding to z = 1) and for a redshifted total mass of M z = 500AfQ. 



1641130]). This evidence includes observed fluxes of radiation that, if isotropic 
and sub-Eddington, would imply IMBH masses above fOO M in many cases; 
thermal peaks at a few tenths of a keV, which, with normal scalings, would 
imply masses of tens to thousands of solar masses; and variability features 
that, if identified with corresponding features observed in systems containing 
stellar-mass black holes, would also suggest IMBHs weighing in at hundreds 
of solar masses. Depending on the threshold one adopts for IMBH likelihood, 
candidates exist in one out of every tens to hundreds of galaxies. Of course, 
the actual number of IMBHs is likely to be far greater than this, in the same 
way that, although we have only about 25 strong candidates for stellar-mass 
black holes in the Galaxy, the total number is likely to be ~ 10 8 . The lack of 
definitive dynamical evidence for IMBHs means, however, that their existence 
is still inconclusive and alternate explanations for the observations have been 
proposed (see [571117] for recent discussions). 

One channel leading to IMBH formation is via the runaway collisions of 
massive stars on time scales too short to allow for stellar evolution, < 3 Myr 
PT] . Recent simulations of runaway collisions with the inclusion of stellar winds 
suggest that winds will prevent the growth of IMBHs in all but the most 
metal-poor environments [41], although these simulations extrapolate wind 
rates from much less massive stars and the collision rates are so high that the 
collision products are likely to be extended bags of gas rather than relaxed 
stars. Alternatively, IMBHs could form through mergers of stellar-mass black 
holes in dense subclusters at the cores of globular clusters; however, recoil 
kicks may eject the products of such collisions from the host globulars [85] . In 
addition to runaway collisions and repeated stellar-mass BH collisions, other 
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IMBH formation channels include IMBH growth through gas accretion early in 
the cluster history [90] and direct collapse of Population III stars (these IMBHs 



could also be seeds for massive black holes, which we discuss in Section 4.1 1. 
We refer the reader to [751174] for reviews. 

Numerical simulations of globular clusters suggest that IMBHs could merge 
with numerous lower-mass compact objects (COs) during the lifetime of the 
cluster PMI771[751[H01[791H41,H51I551 [7D] . through a combination of the emission 
of gravitational radiation, binary exchange processes, and secular evolution of 
hierarchical triple systems. For IMBH mass < 3000M Q , the GWs generated 
during the inspiral of a stellar-mass object (black hole or neutron star, since a 
white dwarf or a main sequence star would be tidally disrupted) into an IMBH 
are potentially detectable by the Einstein Telescope. Ringdown radiation could 
even be detected from more massive IMBHs. 

When the primordial binary fraction in a globular cluster is sufficiently 
high, > 10%, stellar collisions during binary scattering interactions may lead 
to the production of two IMBHs in a single cluster, according to Monte Carlo 
simulations carried out by [46] . Since observations and numerical calculations 
suggest that clusters may be born with large binary fractions (e.g., [55]), the 
formation of two IMBHs may be generic in sufficiently dense and massive 
clusters. If this happens, the two IMBHs will exchange into a common binary 
which, after shrinking via dynamical friction and dynamical encounters with 
other stars, will merge through radiation reaction; all of these processes occur 
on a timescale of < 10 Myr [35] ■ 

Additionally, two IMBHs from different globular clusters may merge during 
the merger of their parent clusters [5J. N-body simulations suggest that, as 
the clusters merge, the IMBHs should form a binary with a peak eccentricity 
of ~ 0.8 [5J, although the residual eccentricity would be negligible by the 
time the frequency reaches ~ 1Hz and the system enters the ET band. The 
binary should merge on a timescale of a few hundred million years through 
a combination of dynamical interactions with stars and gravitational-wave 
emission. If the probability of forming an IMBH in a cluster is high, and if 
clusters merge with a probability ~ 0.1-1 as discussed in [5], the rates of 
mergers of IMBHs originating in different host clusters could be competitive 
with the single cluster channel (see also [TO]). 

In the following two subsections, we estimate the rates with which the 
Einstein Telescope could detect gravitational waves from these two globular 
cluster channels. In Section |3.1| we estimate the rates of detectable inspirals 
of stellar-mass compact-objects into IMBHs (IMRIs) and in Section 3.2 we 
estimate the rate of detectable coalescences of IMBH-IMBH binaries formed 
within a single globular cluster. These subsections contain detailed derivations 
which are presented here for the first time; readers who are not interested in 
the details of these calculations will find the results for IMRIs in Table Q] and 
the results for single-cluster IMBH-IMBH binary mergers in Equation (fl5|. 
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3.1 Intermediate-mass-ratio inspirals into IMBHs 

In an earlier work, a subset of the authors analyzed the possibility of detecting 
intermediate-mass-ratio inspirals (IMRIs) of compact objects into IMBHs with 
Advanced LIGO [70]. In that case, it was found that binary tightening via 
three-body interactions was the dominant channel that led to IMRIs. The 
IMBH, as the most massive object in the cluster, readily switches into compact- 
object (CO) binaries. Once a sufficiently hard CO-IMBH binary is formed, 
the binary will be hardened rather than disrupted by three-body interactions 
with other stars in the cluster. Eventually, as the interacting stars take away 
energy from the binary, the binary will tighten to the point where radiation 
reaction from gravitational-wave emission becomes dominant and drives the 
binary to coalesce. For COs that are neutron-stars or black-holes, it is possible 
to compute the distance to which the gravitational waves can be detected and 
convert this into an estimated detection rate. 

Here, we repeat that calculation for ET sources, with the following two 
major changes. First, we take advantage of the recent development of hybrid 
waveforms that describe all three phases of the coalescence - inspiral, merger, 
and ringdown - to compute the SNR from the full GW signal, rather than just 
the inspiral portion. We use the EOBNR waveforms [23] for this calculation. As 
discussed in Section [2~3] these waveforms have not been tested for mass ratios 
below ^1:4, but, unlike the NSphenom waveforms, they do appear to behave 
correctly in the extreme-mass-ratio limit, ?y — > 0. We use these waveforms as 
there are no better IMRI models available at the present time, but emphasise 
that there will be some corresponding uncertainty in the results, at the level of 
tens of percent. Second, because ET has a low frequency cutoff (~ 1 Hz) that 
is lower than Advanced LIGO (~ 10 Hz), we consider inspirals into 1000 M 
IMBHs along with inspirals into IOOMq IMBHs. We note, however, that for 
higher IMBH masses the IMBH could dominate the dynamics in the center of 
the cluster and a cusp could be formed around the IMBH, possibly increasing 
the importance of the direct-capture scenario [53] ; additional discussion of this 
possibility can be found in Section 2.3 of [70] , 

Approximately (27t/22)M/to* close interactions with stars of mass to* are 
required to reduce the semimajor axis of a CO-IMBH binary with IMBH mass 
M by one e-folding |95| . Stars come within a distance equal to the semimajor 
axis separation, a, of the binary at a rate 



n a M 10 km/s _ x 

10 5.5 pc -3 ioi3 cm 100M© a yr ' 

(2) 

where the bracketed expression is the gravitationally focused cross-section, a 
is the velocity dispersion, and n is the number density of stars in a globular 
cluster, with fiducial values for core-collapsed globulars taken from 94J. The 
last e-folding time dominates the hardening rate, so the hardening time-scale 
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Meanwhile, the gravitational-wave merger timescale for a binary of semi- 
major axis a, eccentricity e, reduced mass approximately equal to the CO mass 
/j, « m, and total mass ps M is [55] 

Tgw ^ 10 m^Iio^J (1 ~ e) yr ^ 10 ^T v T— J 

(4) 

where in the last equality we set e s» 0.98 as the eccentricity after the fi- 
nal three-body encounter, following |45j . Minimizing the total merger time 
Emerge = ^harden + ?gw over a , while setting n, a and to* to their fiducial val- 
ues, allows us to compute the CO-IMBH coalescence rate per globular cluster, 

1 /-^mergc ■ 

To compute the volume within which the Einstein Telescope can detect 



such IMRIs, we follow the procedure outlined in Section 2.2 We use EOBNR 
waveforms and ignore the spin of the IMBH, which we expect to be small, 
S/M 2 < 0.3, after a significant number of minor mergers [69]. We compute 
the range for a "single ET" configuration. The range is a function of the 
rcdshiftcd masses of the IMBH, M z = M(l + z), and the compact object, 
m z = m(l + z). After computing the range, we convert it into a redshift, z, 
by inverting the following expression for the luminosity distance [531 : 



D L (z) = D H (1 + z) I f — 

[Jo [Q M {l + z'f 



(5) 



Here, we implicitly assume a flat universe {fi^ — 0) and use Dm — 0.27, 
f2 A = 0.73, H = 72 km s" 1 Mpc" 1 ; and D H = c/H « 4170 Mpc. We assume 
that the typical source is located near the redshift z that corresponds to the 
search range, and obtain the source-frame masses by dividing the rcdshiftcd 
masses by 1 + z; we use these source- frame masses to compute the merger 
timescale T morgo from Equations ([3| and Q. 

We additionally assume that 10% of clusters form an IMBH and are suffi- 
ciently dense to be relevant to the rate calculation, and that globular clusters 
have a fixed comoving space density of 8.4/i 3 Mpc -3 92 . For h 0.72, this 
yields a density of ~ 0.3Mpc~ 3 for relevant clusters. We compute the comov- 
ing volume up to redshift z by integrating the following expression for dV c /dz 
[S3], with the cosmological parameters defined above: 

=^D* H [n M {i + z? + n A ]- l/2 { f — — d 4——rn) • (6) 



d 



z 



yr, 



The rate of detectable events can then be estimated as ~ 0.3(V r c /Mpc 3 )/[T mcrgc (l+ 
z)], where the factor of (1 + z)^ 1 is included to convert the coalescence time 
measured in the source frame to time measured in the observer frame. 
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100 
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2.2 


31 


3.1 


4 x 10 8 


7 x 10 11 


175 


100 


2 




6.2 


1.0 


51 


1.0 


4 x 10 8 


1.3 x 10 11 


55 


1000 


10 




2.2 


0.4 


710 


7.1 


9 x 10 7 


1.7 x 10 10 


10 


1000 


2 




0.7 


0.15 


870 


1.7 


1 x 10 s 


1 x 10 9 


2 



Table 1 "Single ET" average range, corresponding redshift, source-frame masses, merger 
timescale, comoving volume within range, and detectable event rate for several combinations 
of plausible redshifted CO and IMBH masses. 



Table [T] summarizes the rate predictions for four combinations of M z and 
m z . Although the lack of knowledge about IMBHs and their mass distributions 
makes it impossible to generate firm predictions, and even a lower limit of zero 
IMRIs is possible, it appears that ET may detect hundreds of compact-object 
IMRIs into IMBHs over three years of operation. If ET is operated in the 
xylophone configuration these rates would increase further. 

3.2 IMBH-IMBH inspirals 

In this section, we wish to estimate the rate at which the single cluster channel 
generates IMBH-IMBH binaries that are detectable with the Einstein Tele- 
scope. To do this, we follow the event rate calculation for LISA and Advanced 
LIGO described in [35 . Once a pair of IMBHs is formed in a single cluster, 
they sink rapidly to the center where they form a binary and merge via three- 
body interactions with the stars in the cluster (see [35, TJ for more details). 
Therefore, the rate of IMBH binary mergers is just the rate at which pairs of 
IMBHs form in clusters. The rate of detectable coalescences is 

R _ dN cvcnt = f Mtot '™ x dMtt f\ r m " (Mto " 9) dz d 4 N cvcnt die dVc 
dt J Mtat min t0t Jo Jo dM tot dqdt e dV c dt a dz ' 

(7) 

Here t Q is the time measured in our observer's frame and t e is the time mea- 
sured at the redshift z of the merger; M to t is the total mass of the coalesc- 
ing IMBH-IMBH binary and q < 1 is the mass ratio between the IMBHs; 
Zmax(-/Wtot , q) is the maximum redshift to which the ET could detect a merger 
between two IMBHs of total mass M to t and mass ratio q; dt e /dt a = (1 + z)^ 1 
is the relation between local time and our observed time, and dV c /dz is the 
change of comoving volume with redshift, given by Eq. 
We make the following assumptions: 

— IMBH pairs form in a fraction g of all globular clusters. 

— We neglect the delay between cluster formation and IMBH coalescence, 
since it is expected to be no more than a few tens of millions of years [3 5) . 

— When an IMBH pair forms in a cluster, its total mass is a fixed fraction of 
the cluster mass, M to t — 2 x 10~ 3 M c \. This assumption is based on what 
is typically seen in simulations [47] . As there are no current constraints 
on the mass ratio, we take it to be uniformly distributed between and 
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1. We restrict our attention to systems with a total mass in the IMBH 
range, which we adopt to be between Mtot,iniu = 100M Q and M to ^ max = 
20000M Q . This means we confine our attention to clusters with masses 
5 x 10 4 < M c \/Mq < 10 7 (note that the lower limit is different from that 
chosen in [35] since here we set M to t,min = 1OOM for IMBH sources). 
Thus, 

dM tot dqdt e dV c 9 dM cl dt e dV c 2 x lO" 3 ' () 

The distribution of cluster masses scales as (dN c i/dM c i) oc M~ 2 indepen- 
dently of redshift (122] and the total mass formed in all clusters in this 
mass range at a given redshift is a redshift-independent fraction g c \ of the 
total star formation rate per comoving volume: 

d 2 M SF [ M — d 3 N cl 
Qci = / M c \dM c \, (9) 

Jci dv c dt e y Mclmin dM cl dt e dv c ch 1 ; 

which provides the normalization for dN c \/dM c \: 

d 3 N cl g cl d 2 M SF 1 



(10) 



dM cl dt e dV c ln(M c i, max /M c i, min ) dV c dt e M\ 
— The star formation rate as a function of redshift z is 

d 2 M SF e 3Az \Q M {l + zf + Q A ] l/2 , , , % 

This is the formula used by |107j . in which the star formation rate rises 
rapidly with increasing z to z ~ 2, after which it remains roughly constant. 
As in Section 3.1 we assume a flat universe (S2k = 0), and use flu = 0.27, 
Q A = 0.73, and H = 72 km s" 1 Mpc" 1 . 

Under these asssumptions, Eq. ^ predicts the rate of detectable coales- 
cences per year as 

fl= /n 2 r Xl0 ": 9 / 9cl , / I dq (12) 



ln(Af to t,max/ Mtot,min) 




4M M ,q) e 3Az AitD ^ [ f z dz , 

dz 0.17- 



2 



o [n M {i + z'f + Q A ] 1/2 



Note that here M to t is measured in solar masses and Dh is measured in Mpc. 

Rather than computing z max (M tot ,q) for all values of M tot and q, we rely 
on the following fitting formula for the average range D as a function of the 
redshifted total mass M z = M to t(l+^), obtained by using EOBNR waveforms 
to model the coalescence (see Section 2.2): 

/ , r \ 3 / S 

( #M if M z < M 

D(M z ) = (AMpc){ ) M °{ 11/10 , s-x/2 ^ , (13) 
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where A = 500, M Q = 6OOM for q = 1 and A = 281, M = 450M© for 
g = 0.25. 

Lensing of gravitational wave sources adds some uncertainty to this pic- 
ture. Individual sources can be magnified or de-magnified by lensing, making 
them visible at a greater or lesser distance than predicted by the preceding for- 
mula. Flux conservation ensures that the expected magnification of a source is 
1, meaning no net change in flux (this argument was first elucidated in |117] ). 
which leads us to expect that the total change in the event rate will be small. 
The magnification distribution that arises from lensing peaks at less than 1, 
i.e., a demagnification, and shows an exponential fall-off for large magnifi- 
cations [87 llll6j . The peak moves toward greater demagnification for sources 
at higher redshifts, but this is compensated by a longer tail toward very high 
magnifications. The amount of volume added to a flux-limited sample by these 
highly magnified lines of sight compensates for the smaller volume lost by each 
of the (greater number) of demagnified lines of sight. In [87 , it was shown that 
weak lensing did not significantly change the number counts in a flux limited 
radio sample, but the number of events tended to be increased by the strong 
lensing tail. The total change in number counts was only a few percent. We 
can therefore ignore the effect of weak lensing on the number counts of grav- 
itational wave sources, although it will have an impact on the precision to 
which distances can be measured with gravitational wave observations. 

Ignoring lensing, we can compute z{D\ J ) by inverting Eq. (f5j). For a given 
choice of -M tot and q, the maximum detectable redshift z max (M tot ,q) is then 
obtained by finding a self-consistent solution of 

z(r> h (M to t(l + = Zmax- (14) 



The integrals over M to t and z in Eq. ( 12 ) were evaluated for two specific 
values of q. For q = 1, the total rate was found to be R = 7.5 x 10 4 g g c \; for 
q = 0.25, it was R = 2.7 x 10 4 g g c \. The range varies smoothly with q and so 
we can estimate the integral over q to be the average of these two rates. This 
yields a final estimate of the total rate as 



i? w 500 



(&)(&)*-■ < 15 » 



where we arbitrarily chose g = 0.1 and g c \ = 0.1 as the reference values for 
these unknown parameters. 



4 Sources in low-mass galaxies 

4.1 Light seeds of MBHs at high redshifts 

Supermassive black holes (SMBHs) weighing millions to billions of solar masses 
are nowadays believed to reside in most local galaxies ( |31) and references 
therein). The masses of today's SMBHs exhibit clear correlations with the 
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properties of their host galaxies (luminosity, mass, and stellar velocity dis- 
persion), suggesting there is a single mechanism for assembling SMBHs and 
forming galaxies. The evidence therefore favours a co-evolution between galax- 
ies and SMBHs. 

In the currently favoured cold dark matter cosmology, galaxies today are 
expected to have been built up, via a series of mergers, from small-mass build- 
ing blocks that condensed out at early cosmic times. A single big galaxy can 
be traced back to hundreds of smaller components with individual masses as 
low as ~ 10 5 Mq. Similarly, we expect the SMBHs found in galaxies today to 
have grown partially by accretion and partially by mergers following mergers 
between galaxies (e.g., |112| l68]). so that a single SMBH can be traced back to 
some number of 'seed' black holes at early times [113] . There are large uncer- 
tainties in this picture, however. Did seed black holes form efficiently in small 
galaxies (with shallow potential wells) at early times, or was their formation 
delayed until substantial galaxies with deeper potential wells had been formed? 
This is a key question, as the mass and the occupation number of the seeds 
ultimately dictates the occupation number of SMBHs in galactic centers. 

The formation of SMBHs is far less well understood than that of their light 
stellar-mass counterparts. The 'flow chart' presented by [M] still stands as a 
guideline for the possible paths leading to the formation of SMBH seeds in 
galactic nuclei. One possibility is that the seeds of SMBHs were the remnants 
of the first generation of stars, formed out of zero-metallicity gas 65j. In a cold 
dark matter universe, structure builds up hierarchically, so the smaller clumps 
at the earliest cosmic times have shallower potential wells. Stars cannot form 
until the clumps are sufficiently big to provide a potential well deep enough 
to pull in gas that can cool radiatively and contract to make a protostar. 
This requires dark matter clumps - minihalos - of ~ 10 6 Mq at redshifts of 
z ~ 20. The first stars forming in these minihalos develop under very differ- 
ent conditions from present-day stars: there are no heavy elements (so that 
molecular hydrogen is the only effective coolant), no dust, and no magnetic 
fields. These conditions mean that these 'Population III' stars were likely very 
massive, having characteristic masses of the order of ~ 100 M Q (e.g., [2T3I 
I5T1I21H20 )■ This prediction relies on the absence of efficient cooling agents in 
the primordial metal-free gas. If Population III stars form with masses 40 
M < M <140 M Q or M> 260 Mq, they are predicted to collapse and form 
IMBHs directly with little mass loss [37], i.e., leaving behind seed IMBHs with 
masses Mbh ~ 10 2 — 10 3 Mq. This is a plausible formation mechanism for 
the seeds upon which supermassive black holes are grown |112j , although more 
massive black holes may have been formed after the epoch of the first stars in 
dark-matter halos with virial temperatures of ~ 10 4 K 21 . K)(i. 1"> via 'direct 
collapse', as described in the following. 

Direct collapse models for MBH formation rely on the collapse of supermas- 
sive objects formed directly out of dense gas [4"%, 63,22, 60l lTM6"2"] . The phys- 
ical conditions (density, gas content) in the inner regions of mainly gaseous 
proto-galaxies make these loci natural candidates, because the very first proto- 
galaxies were, by definition, metal-free, or at the very least very metal-poor. 
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Enriched halos have a more efficient cooling, which in turn favours fragmenta- 
tion and star formation over the efficient collection of gas conducive to MBH 
formation. In a typical galaxy, however, the tidally induced angular momen- 
tum would still be enough to provide centrifugal support at a distance ~ 20 
pc from the centre, and halt collapse, ultimately leading to the formation of 
a disk. Additional mechanisms inducing transport of angular momentum are 
needed to further condense the gas until conditions fostering MBH formation 
are achieved. An appealing route to efficient angular momentum shedding is 
by global dynamical instabilities, such as the "bars-within-bars" mechanism, 
that relies on global gravitational instability and dynamical infall [1041115) . 
Self-gravitating gas clouds become bar-unstable when the level of rotational 
support surpasses a certain threshold. A bar can transport angular momen- 
tum outward on a dynamical timescale via gravitational and hydrodynamical 
torques, allowing the radius to shrink. Provided that the gas is able to cool, 
this shrinkage leads to even greater instability, on shorter timescales, and the 
process cascades. This mechanism is a very attractive candidate for collecting 
gas in the centres of halos, because it works on a dynamical time and can 
operate over many decades of radius. It has also been proposed that gas ac- 
cumulation in the central regions of protogalaxies can be described by local, 
rather than global, instabilities. During the assembly of a galaxy disc, the disc 
can become self-gravitating. As soon as the disc becomes massive enough to 
be marginally stable, it will develop structures that will redistribute angular 
momentum and mass through the disc, preventing the surface density from 
becoming too large and the disc from becoming too unstable. To evaluate the 
stability of the disc, the Toomre stability parameter formalism can be used 
[501162) . The gas made available in the central compact region can then form 
a central massive object, for instance via the intermediate stage of a 'super- 
massive' star [55,13 , or a 'quasistar', an initially low-mass black hole rapidly 
accreting within a massive, radiation-pressure-supported envelope, [15, 3]. In 
both cases, the mass function of seeds is predicted to peak at 10 5 — 10 6 A/ Q 

As described in Section 3.1 the formation of an IMBH as a result of dynam- 
ical interactions in dense stellar systems is a long-standing idea, which could 
also create intermediate mass MBH seeds. This process could have been very 
effective in the very first stellar clusters that formed in high-redshift proto- 
galaxies, when the Universe was not as metal-rich as now. Low metallicity 
favors the growth of a very massive star, the precursor of an IMBH remnant. 
The mass loss due to winds is significantly reduced in metal-poor stars, which 
greatly helps in increasing the mass of the final IMBH remnant (cf. |121j ). The 
formation of stellar clusters and the possible evolution of the stellar systems up 
to IMBH formation are explored in [55]. Figure [3] shows three mass functions 
for three different MBH 'seed' scenarios: direct collapse [TS], runaway stellar 
mergers in high-redshift clusters, and Population III remnants [112] , 

It is uncertain how many MBH 'seeds' formed, and in which mass range. 
Equally uncertain is how these 'seed' black holes grew within their host mini- 
halos. It is not obvious if efficient accretion onto these seeds could have taken 
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Fig. 3 Mass function of seed MBHs for three different formation scenarios. Left: direct 
collapse I113| : centre: runaway stellar mergers in high-redshift clusters 28 ; right: Population 
III remnants 1651 . Note the different y-axis scale for the Population III case. 



place, at least early on, in the fragile environment that the shallow potential 
wells of minihalos represent [75]. It is likely that seed IMBHs can grow ef- 
ficiently only if they are hosted in the most massive galaxies at these early 
cosmic epochs, while IMBHs in an 'average' galaxy could have experienced 
intermittent and inefficient accretion, thus leaving behind a population of 
underfed IMBHs with a mass range similar to that of the original seeds, 
M BH ~ 10 2 - 1O 3 M . 

The Einstein Telescope will be able to probe mergers between black-hole 
seeds at high redshift, and thus help to distinguish between these various 
channels for seed formation. An estimate for the ET event rate under the 
Pop III model can be computed using Monte-Carlo merger-tree realizations 
based on the extended Press-Schechter formalism [S3], as described in [1121 
1114] . This was done in [100], and the results which we now quote are all taken 
from that work and the companion paper |39| . These papers considered four 
different models that were based on the same merger tree realisations (taken 
from |112U114| ). but differed in the initial mass distribution of seeds and in 
the prescription for accretion onto the seed black holes. In these scenarios, 
which were all based on having the light, Pop III, remnants described earlier 
in this section as the seeds for black hole formation, a single ET would detect 
~ 1-10 seed mergers, depending on the model. The detected mergers would 
be between black holes with total mass ranging from 2Af m i n up to ~ 1000 M@, 
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where M m ; n is the mass of the lightest seed black hole in the initial mass 
distribution. This minimum seed mass is rather uncertain and depends on 
the details of the model used, as discussed earlier. In the scenarios considered 
in |100j . M m j n was either 1OM or 15OAf . The detected events would be seen 
at redshifts z ~ 1-7, although this could extend to z ~ 12 for the lightest 
seed model, which had M min = 1OM . If ET was operated in the xylophone 
configuration described in Section |2.1| the number of events seen would be 
increased to several tens, and these would be out to a redshift z ~ 15 [39] , 

Figure [4j reproduced from the data in |39| . shows how the number of events 
seen by ET over three years varies as a function of the signal-to-noise ratio 
required in a single 10km right-angle interferometer for detection. The SNR 
required in the network of detectors is likely to be ~ 8, although this depends 
somewhat on data-analysis issues, and on the amount of source confusion in 
the data stream. A network SNR of 8 corresponds to an SNR in the single 
right-angle detector of 5.3 for a single ET, or SNRs of 4.8, 3.9, 3.8 and 3.1 
for the network configurations (i) - (iv) described in Section 2.1 In Figure [4] 
we show results for two of the four light-seed models considered in 100 , and 
for both the baseline and xylophone configurations of the detector. The rate 
for the baseline ET configuration is rather sensitive to the SNR that is ul- 
timately required for a confident detection, but the xylophone configuration 
is more robust, as it has improved sensitivity at just the right frequency for 
systems with mass in the 1OO-1OOOM range. The mergers seen by ET will be 
complementary to mergers between heavier black holes that will be seen by 
space-based detectors such as LISA, ALIA or DECIGO [39]. The combination 
of detectors will provide a nearly complete survey of mergers between galactic 
black holes, yielding important constraints on astrophysical models of galaxy 
formation and growth. The utility of observations of binary black hole sys- 
tems with multiple detectors was also discussed in [10] , with specific reference 
to IMBH-IMBH binary mergers arising from the mergers of globular clusters 
containing central IMBHs. They found that, if LISA and ET were operating 
concurrently, the same IMBH binary could be detected by both detectors with 
a time separation of a few months. Even if observations of the same system are 
not made with different detectors, each detector will provide a measurement 
of the rate of black hole mergers in a different black hole mass range. These 
rate observations over the black hole mass spectrum will provide important 
constraints on models of black hole growth. 

One important question is whether ET will be able to distinguish be- 
tween black-hole mergers coming from this channel, and those described in 



Section 3.2 that arise in globular clusters. To provide constraints on merger 



histories, it is necessary to know that an observed event is associated with a 
galaxy merger. The masses and redshifts of events will provide some informa- 
tion, but more work is required to understand what observational signatures 
provide the best discriminating power. We would expect mergers between seed 
black holes to occur over a range of redshifts, with some events at redshifts 



> 10. In the mechanism described in Section 3.2 the black-hole binaries 



form and merge very quickly, so this could also produce events over a range 
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Fig. 4 Number of events detected by the Einstein Telescope in three years, as a function 
of the required signal-to-noise ratio threshold in a single right-angle detector. Results are 
shown for both the baseline and xylophone configurations of the Einstein Telescope, and 
for two different astrophysical models — Volonteri-Haardt-Madau (VHM) with equal mass 
seeds (VHM, ems) and VHM with a seed mass distribution (VHM,smd). Details on these 
models can be found in 1391 . 



of redshifts. However, the distinction between these two formation channels 
becomes increasingly vague at high redshift, when galaxies are in the process 
of formation. What is important for the light-seed scenario is that black holes 
of low mass, ~ IOOMq, exist at high redshift. Therefore, being able to identify 
an event as being between two ~ 100M Q black holes at redshift z > 5 would 
be an important constraint, regardless of how that seed had initially formed. 
A single ET cannot measure the six extrinsic parameters of a merger source 
on its own — at least one additional non-colocated detector will be required. 
Possible network configurations were discussed in Section [iOj With one addi- 
tional 10km detector at the location of LIGO Hanford, the ET network will be 
able to determine the luminosity distance of a source to an accuracy of ~ 40%. 
Adding a third 10km detector at the site of LIGO Livingston or upgrading the 
detectors to ETs improves this modestly to ~ 30% 100,39 . There will also be 
an additional distance error due to weak lensing of the signal, but this will be 
considerably smaller (~ 10%, see, for example, |103] ) than the intrinsic error 
from the gravitational wave observations. If we assume that the luminosity 
distance is converted into a redshift using the concordance cosmology at that 
time, the redshift error will be comparable to the distance error. Thus, an ET 
network should be able to say with confidence if an event is indeed occurring 
between two ~ 100M Q black holes at high redshift, z ~ 5. 
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4.2 MBHs in dwarf galaxies 

There are two simple arguments that lead us to believe that ~ 10 2 -10 3 Af Q 
black holes might inhabit the nuclei of dwarf galaxies today. Firstly, the mass 
of SMBHs detected in neighboring galaxies scales with the bulge mass — or 
stellar velocity dispersion (A/bh — c) — of their host galaxy II) 32[I~)0I. 
The lowest-mass galaxies currently known have velocity dispersions a ~ 10- 
20 km s _1 |115j . If we extrapolate the Mbh — cr correlation to these a values, 
we expect the putative IMBHs to have masses in the range of hundreds to 
thousands of solar masses. 

Secondly, as SMBHs grow from lower-mass seeds, it is natural to expect 
that a leftover population of progenitor IMBHs should also exist in the present 
universe. Indeed, one of the best diagnostics of 'seed' formation mechanisms 
would be to measure the masses of IMBHs in dwarf galaxies. This can be un- 
derstood in terms of the cosmological bias. The progenitors of massive galaxies 
have a high probability that the central SMBH is not "pristine", that is, it 
has increased its mass by accretion, or it has experienced mergers and dynam- 
ical interactions. Any dependence of A/bh on the initial seed mass is largely 
erased. However, low-mass galaxies undergo a quieter merger history, and as 
a result, at low masses the BH occupation fraction and the distribution of BH 
masses still retain some "memory" of the original seed mass distribution. The 
signature of the efficiency of the formation of SMBH seeds will consequently 
be stronger in isolated dwarf galaxies [113] . 

One hopes that the next generation of 25-30m optical/IR telescopes operat- 
ing at their diffraction limit (~ 4 milliarcsec) can provide the first constraints 
on the presence of IMBHs in dwarf galaxies, but the detection of gravitational 
waves from a central IMBH in a dwarf galaxy undergoing a merger is possibly a 
more promising probe. Dwarf galaxies have a very quiet merger history, hence 
we do not expect many IMBH-IMBH mergers involving dwarf galaxies at the 
present epoch, or in the low-redshift universe. The seed black hole mergers 



discussed in Section |41] probe a separate population of mergers, between the 
progenitors of galaxies which are more massive today. However, gravitational 
waves may also be generated in dwarf galaxies by mergers between the central 
IMBH and stellar remnants in the centre of the dwarf. These are analogous to 



the globular-cluster IMRI sources described in Section [31 

We can derive an estimate of the event rate based on the expected num- 
ber of dwarf galaxies which can possibly host IMBHs in the interesting mass 
range. Theoretical models of SMBH formation and evolution, where the seeds 
of MBHs arc Population III remnants |112j . can be used to look for the dis- 
tribution of IMBHs in dwarf galaxies. Using the dynamical model of [110) . we 
estimate a number density of IMBHs, t^imbh ~ 0.02-0.1 Mpc -3 . 

When we calculate the event rate of BH-IMBH mergers in dwarf galaxies, 
we have to further correct for the fact that only a small fraction of these tiny 
satellites do indeed form stars [H]. Based on 42], we estimate that a fraction 
/* = 0.1 — 0.2 of dwarfs in the w sa t ~ 10 — 20 km s _1 range formed stars 
(which will eventually leave behind stellar mass BHs that can merge with 
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the central IMBH). The number density of IMBHs that can be ET sources is 
therefore uet = f* "imbh ~ 0. 001-0. 06Mpc -3 . This number density is about 
an order of magnitude lower than the number density of globular clusters used 
to normalise the rates in Section 13.11 

The capture mechanisms that lead to IMRIs in dwarf galaxies are likely to 
be the same as those that operate in globular clusters. The event rate for the 
binary-hardening mechanism scales with the stellar density, n, as n 4 / 5 , while 
the other mechanisms, such as direct capture, should scale approximately with 
n (this is the same n that enters Eqs. (|2])-(|3|). The core stellar densities in 
nearby dwarf galaxies are typically much lower than in core-collapsed globular 
clusters, e.g., the estimate for Fornax is ~ 10 _1 pc -3 [73] and for Sagittarius 
is ~ 10- 3 pc- 3 [67], compared to - 10 5 - 5 pc" 3 for globulars [94]. The IMRI 
rates for dwarf galaxies are thus likely to be orders of magnitude lower than 
those for globular clusters. Therefore, although it is not inconceivable that 
ET will detect events from dwarf galaxies, any events would be serendipitous. 
Moreover, there are no obvious characteristics which would allow an observer 
to distinguish between an event in a dwarf galaxy from one in a globular 
cluster based on the GW signature alone. Nonetheless, it might be possible 
to make qualitative statements about dwarf galaxy IMBH populations. For 
instance, if ET does not detect any mergers between seed black holes at high 



redshift, of the type described in Section 4.1 it is very likely that BH seeds 



were heavy and not light. This would suggest dwarf galaxies would not contain 
light leftover BH seeds, and consequently that all of the observed IMRIs are 
occurring in globular clusters. Similarly, if seed mergers are detected but the 
rate of IMRIs is low or zero, it might suggest that IMBH formation in globular 
clusters is inefficient and any observed IMRIs are in dwarf galaxies. More 
refined modelling and detailed calculations are needed to understand/prove 
the robustness of these expectations, especially in view of the small number of 
seed black hole merger events and dwarf galaxy IMRIs that are predicted. In 
summary, while the dwarf galaxy channel should not be ignored completely, 
it is very unlikely to be a significant contributor to ET events or science. 



5 Speculative sources 

In this section we discuss some more speculative sources that might be ob- 
served by a future low-frequency ground-based interferometer such as ET. We 
examine first the possibility of observing orbiting or rotating white dwarfs near 
the high end of allowed masses, then discuss how the observation of eccentric 
compact binaries could illuminate their dynamical origin. 



5.1 Orbiting white dwarfs 

A gravitationally bound object of average density p has a maximum orbital, 
rotational, or acoustic frequency / max oc (G7?) 1 / 2 . For neutron stars this max- 
imum is ~ 10 3 Hz. White dwarfs are much more extended objects, but near 
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their maximum masses their densities are sufficient to reach / max ~ 1 Hz. For 
example, from the classic work [49], a magnesium white dwarf with maximum 
mass M max = 1.363 M has a radius R = 2.57 x 10~ 3 R = 1.79 x 10 s cm and 
therefore (Gp) 1 ' 2 = 2.7 Hz. In the few-Hz range, therefore, one will potentially 
see gravitational waves from the most massive white dwarfs. 

If we consider specifically such a white dwarf in a binary orbit with a 
neutron star, black hole, or another white dwarf, then the orbital frequency 
at the point of tidal disruption of the dwarf depends weakly on the mass of 
the companion. For example, suppose that the equilibrium mass and radius 
of the white dwarf are respectively Mwd and i?wD, and that the companion 
is a compact object of mass M comp . When the orbital separation a is a ~ 
2i?wD(A^compA%D) 1 ' 3 , tidal stripping begins [119| and others]. The orbital 
frequency at this point is 

w = ^G(M comp + M WD )/a 3 ~ 0.7(1 + Mw D /Mcomp) 1/2 (Gp) 1 / 2 . (16) 

The gravitational wave frequency is /gw = 2/ or b = implying a maximum 
frequency of ~ 1 Hz for comparable-mass objects such as a neutron star and 
a heavy white dwarf, and a maximum that is ~70% of this if the companion 
is a much more massive object such as an IMBH. 

We have relatively few candidates for massive white dwarfs, hence al- 
though there is a significant literature related to lower- frequency radiation 
from white dwarf binaries (e.g., [29j and many subsequent papers) their num- 
bers are difficult to estimate (see |lllj for a recent discussion). Models of 
the mass distribution suggest that perhaps ~ 0.1 — 1% of white dwarfs have 
masses near Mwd = M max (e.g., see figure 10 of 26 ). Our requirement that 
both white dwarfs have masses near the maximum means that the mass ratio 
is greater than 2/3, and thus there will be a merger instead of stable mass 
transfer (see [75]). If we estimate that that there are 2.5 x 10 8 double white 
dwarf systems in a galaxy like the Milky Way [82], and that ~ 50% of the 
massive ones have semimajor axes that allow merger by gravitational radia- 
tion within 10 10 yr (corresponding to the ~ 48% merger fraction from [82"]). 
then we expect massive white-dwarf binaries to merge at a rate per galaxy of 
- (0.001 - 0.01) 2 x 0.5 x 2.5 x 10 8 /10 10 yr" 1 - 2 x 10~ 9 - 10~ 6 yr" 1 . At the 
high end this is similar to the low end of NS-NS merger rate estimates [S3] • If 
the ET is sensitive to such mergers out to ~ 200 Mpc, which may be optimistic 
given their low GW frequencies, one event per few years could be detected. 
Detection of these events would indicate rather precisely the maximum aver- 
age density of white dwarfs, and would thus be a mechanism for establishing 
their mass-radius relation near the maximum mass. 



5.2 Rotating hypermassive white dwarfs 

Another possibility, suggested to us by [86 , is that two white dwarfs with 
more typical masses Mwd < 1 M might merge in a binary and produce a 
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hypermassive white dwarf that spins rapidly enough that it is deformed into 
an ellipsoid. This is a promising candidate to explain some fraction of Type la 
supernovae [9"8] . 

To evaluate this prospect we note that if a Newtonian perfect fluid (a good 
model for a white dwarf) rotates uniformly then above a certain critical angular 
momentum, L cl n, for a given mass, M, the equilibrium configuration splits 
off from the axisymmetric Maclaurin spheroids (which emit no gravitational 
radiation) to the Jacobi ellipsoids. If the three axes of the ellipsoids are a 3 < 
&2 < a>i, then, according to [37J section 39], the critical angular momentum is 

L crit « 0.3(GAf 3 a) 1/2 (17) 

where a = (a^as) 1 / 3 . If two white dwarfs both of mass M/2 and radius R 
spiral slowly together, then their angular momentum at the point of contact is 
L = n-sJlGMR = \/2/4(GAf 3 i?) 1 / 2 = 0.35(GM 3 i?) 1 / 2 . Since the equilibrium 
radius of the hypermassive object is smaller than the radii of the original white 
dwarfs, the angular momentum is sufficient to produce an ellipsoidal figure. 
Again from j2jj section 39], the angular velocity of this configuration will be 
fl f=a (Gp) 1 / 2 and hence the dominant gravitational wave frequency will be 
/gw « (Gp) 1/2 A- 

The amplitude of gravitational waves depends on the ellipticity e={I\ — 
I2) / 'I3, where Ii indicates the moment of inertia along axis i. Near the critical 
angular momentum, slight changes in L produce large changes in e, and e of 
several tenths is possible. Gravitational waves remove rotational energy from 
the star, such that 

u, = -f^ e 2 / 3 c 5 (18) 
c 

where uj = tt/gw- As a result, the characteristic spindown time is 

Spindown = | « 200 yr I— I I — jiYUz) ' ^ ' 

The sweep rate at 1 Hz is, therefore, ~ 1 Hz/200 yr ~ 10 -10 Hz s _1 . For 
an integration of ~ 10 5 s the frequency would stay in a single frequency bin 
of Af = 1/10 5 s = 10~ 5 Hz. Since the spindown rate will remain constant 
for a much longer time a search for a simple linear drift may make practical 
integrations over weeks to months. This would partially offset the low expected 
amplitudes. For comparison, continuous wave searches in LIGO are routinely 
done for spindown times as low as ~ 1000 years at frequencies of ~ 100Hz [T], 
so a search for spindown times of ~ 200 years at / ~ 1Hz is certainly feasible. 

Type la supernovae are estimated to occur once per 1000 years in galaxies 
such as the Milky Way [55] > so even if only 1-10% of SNe la are binary mergers, 
the overall astrophysical rate is competitive with double neutron star mergers. 
Even though the ET sensitivity to gravitational waves from these binary white 
dwarf mergers will be much lower than for double neutron star mergers, the 
detection of gravitational waves from any such event may provide a new view 
on these important supernovae. 
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5.3 Eccentric binaries 



In the sensitivity bands of second-generation gravitational wave detectors such 
as Advanced LIGO and Advanced Virgo, most compact binaries will be very 
close to circular. (For a proposed scenario where this may not be true, see 
[S3] : another possibility includes direct captures of compact objects by IMBHs 
as precursors to eccentric IMRIs in globular clusters, although this formation 



mechanism is uncommon relative to the one described in Section 3.1 which will 



produce circular IMRIs). This is because for moderate to high eccentricities, 
gravitational radiation essentially reduces the semimajor axis of a binary while 
keeping the pericenter fixed. Therefore, to have palpable eccentricity at a given 
frequency, the pericenter at formation or at the last dynamical interaction must 
be inside the radius of a circular orbit at that frequency. For example, a binary 
of two 1.4 Mq neutron stars must have a pericenter less than 700 km to be 
significantly eccentric at a gravitational wave frequency /gw = 2/ or b = 10 Hz. 
This is highly improbable for a field binary, and is even difficult to arrange for 
binary-single scattering in dense stellar environments. 



Somewhat higher eccentricities can be obtained via the Kozai secular res- 
onance [ST]. As explored in the context of black holes by [TrUlllS] . a binary- 
binary interaction can result in a stable hierarchical triple in some tens of 
percent of encounters. If the inner binary and the outer tertiary have orbital 
planes that are inclined significantly with respect to each other, then over many 
orbital periods the inclination and eccentricity of the inner binary change pe- 
riodically, leading at points in the cycle to very small pericenters and thus 
potentially observable eccentricity after the gravitational-wave driven inspi- 
ral. The eccentricity at 40 Hz is almost always very small (below 0.1), but at 
10 Hz there are a few orientations in which the eccentricity can be a few tenths 
[118j . At still lower frequencies the eccentricity will be yet higher, because for 
low eccentricities, e, e oc J -19 / 18 . 



The preceding discussion implies that detector sensitivity at low frequencies 
will be important to determine the origin of compact binaries. In-situ forma- 
tion from a massive main-sequence binary is still highly unlikely to produce 
detectable eccentricities: in order to have eccentricity at 1 Hz, the pericenter 
distance would have to be < 3000 km immediately after the second supernova. 
In contrast, dynamical effects such as the Kozai resonance are expected to pro- 
duce eccentric orbits at a few Hz. As a result, observation of a few BH-BH or 
BH-NS inspirals at a few Hz will illuminate their formation processes in a way 
that is not as easy at higher frequencies. We note, however, that simulations 
such as those in [55] suggest that of the few per year to few tens per year of 
black hole mergers in globulars that are expected to be seen with Advanced 
LIGO, less than 10% are initiated by the Kozai process. The greater reach 
of the Einstein telescope will enhance the total numbers, but binaries with 
palpable eccentricity in the ET band are still expected to be a minority. 
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6 Scientific impact of ET observations 

ET detections of any of the systems described in this paper will yield important 
science products, which we now discuss. 

6.1 Astrophysics 

The very existence of BHs in the 100-1000M Q range is uncertain, so a single 
robust detection of an IMBH by ET will be of huge significance. If ET detects 
any seed black hole mergers at high redshift, it will be strong evidence that 
black hole seeds were light, which will help discriminate between light and 
heavy seed scenarios for the growth of structure in the Universe. Observations 
of mergers between more massive black holes with LISA do not have the same 
discriminating power, as they cannot distinguish between ~ 10 5 A/ Q MBHs 
that formed through direct collapse or the collapse of a massive Pop III star 
and those that formed through a sequence of mergers |102[ 39 . A significant 
number of ET detections of seed black hole mergers may provide constraints 
on the mass distribution of black hole seeds, and their early accretion history. 

Detection of a significant number of IMRIs with ET will indicate that 
IMBHs form readily in globular clusters (since the rate of IMRIs in dwarf 
galaxies is so low). The characteristics of the IMRI events will provide con- 
straints on the astrophysics of dense stellar environments, and on the efficiency 
of capture processes operating within them. 

If ET detects white dwarfs undergoing tidal disruption, it will provide 
important constraints on the physics of degenerate matter, including the max- 
imum density and mass that white dwarfs can reach. Detections of rotating 
hypermassive white dwarfs would provide information about proposed chan- 
nels leading to supernovae. Finally, the detection of a significant population of 
eccentric coalescing binaries will shed light on the efficiency of the processes 
that drive eccentricity growth in binaries, such as the Kozai mechanism. 



6.2 Fundamental physics 

ET IMRI sources can also be used for testing aspects of relativity theory, in 
particular verifying that the central object is indeed a black hole as described 
by the Kerr metric of general relativity. This has been explored extensively in 
the context of extreme-mass-ratio inspiral events detectable by LISA (see, for 
example, [9] and references therein) . In the course of an inspiral, the orbit of 
the smaller object traces out the spacetime geometry of the large body and 
hence the emitted gravitational waves encode a map of the spacetime struc- 
ture. One way to characterize this is in terms of the multipole moments of the 
spacetime. It was demonstrated by Ryan [99], for nearly circular and nearly 
equatorial orbits, that successive multipole moments of an arbitrary space- 
time are encoded at different orders in an expansion of the orbital precession 
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frequencies as functions of the orbital frequency. Since these frequencies can 
be measured from the emitted gravitational waves, a multipole map of the 
spacetime can in principle be measured. Similar multipole measurements are 
also possible from observations of ringdown radiation following mergers |18j . 
For a Kerr black hole, the mass, M, and angular momentum, S, determine all 
higher-order mass, Mi, and current, Si, multipole moments of the spacetime: 

Mi+iSi = M(iS/M) 1 . (20) 

Measuring just three multipole moments and finding them to be inconsistent 
with this formula is therefore enough to demonstrate that the central object 
is not a Kerr black hole. 

For IMRIs, it has been shown that Advanced LIGO could measure an 
O(l) fractional deviation in the mass quadrupole moment, M2, for typical 
systems [23]. Corresponding results have not yet been computed for ET. How- 
ever, ET will improve this significantly for two reasons — (i) the SNR of a 
source at fixed distance will increase by a factor of 10 or more; and (ii) ET 
will observe the sources at lower frequencies. The ability to measure multipole 
moments improves significantly with the number of gravitational-wave cycles 
observed. At the leading-order Newtonian approximation, a 1M© + lOOAf© 
system has ~ 500 cycles remaining until plunge when the frequency is 10Hz, 
but this increases to ~ 1500 for a frequency of 5Hz, ~ 4000 for 3Hz and 
~ 25000 for a frequency of 1Hz [33]. ET should thus be able to carry out 
tests of the Kerr nature of the central object that are significantly better than 
those possible with Advanced LIGO. Further research is required to quantify 
the improvement that will be possible, and how this will compare to expected 
results from LISA EMRI events. 



6.3 Uncertainties 

There are various uncertainties which will affect the scientific impact of ET 
measurements discussed above. One important consideration is how to distin- 
guish between IMBH events that arise from seed black holes and those that 
arise from IMBHs formed in globular clusters. Using ET measurements to 
constrain hierarchical structure formation relies on identification of mergers 
as seed black hole mergers, but, as we have seen, there may also be IMBH 
binary mergers in globular clusters. The masses and redshifts of the events 
may provide a robust discriminator, but more work is needed to understand if 
this is indeed the case, or whether other characteristic features exist that can 
be exploited. 

The eventual sensitivity that is achieved by ET also has bearing on these 
results. The speculative sources that were discussed in Section [5] rely on ET 
having sensitivity in the l-10Hz band, and low-frequency sensitivity should 
also improve the accuracy of tests of relativity using IMRIs. ET may only have 
sensitivity down to a frequency of ~ 3Hz, which will impact all of this science 
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and perhaps eliminate the possibility of detecting gravitational radiation from 
massive white dwarfs. This must be properly quantified in the future. 

Finally, there are open questions regarding ET data-analysis. The ET data 
stream will be very source-rich, and so the identification of individual sources of 
different types in the presence of this confusion will be a challenging problem. 
For instance, neutron star binary systems will create a confusion background 
near 1Hz jJT]. The data-analysis challenges for ET will inevitably change the 
SNRs required for detection of individual sources and therefore the rate pre- 
dictions, and the accuracy with which source parameters can be estimated. 
However, the rate uncertainties arising from the data analysis will most likely 
be small compared to the order-of-magnitude uncertainties that are present in 
the astrophysical rate predictions. 



7 Summary 

We have discussed gravitational waves generated by intermediate-mass black 
holes as possible sources for the Einstein Telescope. Intermediate-mass black 
holes may be formed via two alternative channels — (i) they may be formed 
in the early Universe if MBH seeds are light (seed IMBH) ; (ii) they may form 
in globular clusters via runaway collisions between stars (cluster IMBH). In 
both cases, there are two distinct types of system that might be sources of 
gravitational waves for ET — (a) mergers between binaries containing two 
IMBHs; (b) mergers of stellar remnants with IMBHs (IMRIs). 

Mergers between seed IMBHs occur following galaxy mergers during the 
hierarchical assembly of structure. If MBH seeds are light, ET could detect a 
few to a few tens of seed black hole merger events over three years at redshifts 
as high as z ~ 8-10. An ET network would, in addition, be able to determine 
the luminosity distance to these events to an accuracy ~ 30%, which is suffi- 
cient to say confidently that an event involves intermediate-mass black holes 
and is occurring at high redshift. IMRIs involving seed IMBHs could occur in 
dwarf galaxies, but the event rate is probably very low, which makes it un- 
likely that this will be a significant contributor to the ET event rate. If cluster 
IMBHs form readily, binary IMBHs in globular clusters might be detected by 
ET at a rate of ~ 500 per year. Core-collapsed globular clusters are also a 
more promising host for IMRIs detectable by ET and the IMRI event rate for 
ET could be as high as a few hundred per year. However, there are significant 
uncertainties, not least of which is whether IMBHs form at all in the stellar 
environments of globular clusters. 

The improved sensitivity of ET at low frequency may also allow the detec- 
tion of several speculative sources. High-mass white dwarfs can survive tidal 
disruption long enough to reach orbital frequencies / max ~ 1Hz in binaries. 
Hypermassive white dwarfs formed by the mergers of normal white dwarfs in 
binaries could also be sources for gravitational waves at frequencies around 
1Hz as they will be rapidly rotating and can support relatively significant el- 
lipticities. ET could detect these two types of event at a rate of one per few 
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years, but this number is extremely uncertain. Finally, dynamical processes 
such as the Kozai mechanism can excite sufficiently high eccentricities in BH- 
BH and BH-NS binaries, that there would be significant residual eccentricity 
when their orbital frequency is in the l-10Hz range that ET will probe. These 
systems would circularize before reaching orbital frequencies in the Advanced 
LIGO band. ET might detect several eccentric binaries per year, but this rate 
depends on the fraction of binaries with residual eccentricity and as yet un- 
known details of the ET data analysis. 

ET detections of any of these sources would have significant impact on our 
understanding of various astrophysical processes, as well as being useful for 
fundamental physics. 
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